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g I We have used x-ray anomalous diffraction to extract the x-ray structure factor of 

■ InAs quantum stick-like islands, embedded in InP. The average height of the quantum 

O ■ sticks (QSs), as deduced from the width of the structure factor profile is 2.54nm. The 

InAs out of plane deformation, relative to InP, is equal to 6.1%. Diffraction Anoma 
^ ' lous Fine Structure provides a clear evidence of pure InAs QSs. Finite Difference 

Method calculations reproduce well the diffraction data, and give the strain along 
the growth direction. Chemical mixing at interfaces is at most of IML 
PACS : 68.65.La 61.10.Nz 



Band structure gives the optical and electronic properties of materials. It can be strongly 
modified by reducing the size of semiconductor material down to a length scale comparable 
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to the effective wave-lengtli of the carriers, i.e. in the order of several nanometers, leading 
to discrete energy levels Much study has been done on InAs quantum structures such as 
Quantum Wires (QWrs) and Quantum Dots (QD). With an emission wavelength which can 
be around 1.55^m, such nanostructures are of interest for the next generation of integrated 
circuits (telecommunication relevant range of 1.3-1. 6//m). The nano-objects we are interested 
in are grown by Molecular Beam Epitaxy (MBE) and obtained via the Stranski-Krastanov 
growth. The lattice mismatch "^"^"""^"^ is about 3.2%. Recent studies yl have shown 
that a strong stress anisotropy appears during the growth deposition leading to a higher 
stress along the [110] direction than along the [110]. The stress is released first in the [110] 
direction leading to the QWrs formation. To be suitable for devices, the nanostructures 
are encapsulated with InP or embedded in a superlattice, they must be homogeneous in 
size, shape and composition, to provide well defined emission wavelengths. The knowledge 
of strain field, chemical gradients, chemical mixing at the interface, is of great importance 
to understand the growth dynamics as well as the electronic and optical properties of the 
nanostructures. 

We report on an x-ray study of InAs stick-like islands embedded in InP with 10 nm thick 
capping layer. These samples are obtained by optimisation of the MBE growth parameters 
to minimize the As/P exchange and reduce the height dispersion of the InAs islands 3,1^. 
The width of photoluminescence (PL) peaks reveals a very good height homogeneity p[, as 
confirmed by Transmission Electron Microscopy (TEM) measurements. The InAs islands 
obtained with a nominal deposition thickness of 4ML are 50 to 200nni long along the [110] 
direction ([T^). As shown by cross-section TEM image the QSs exhibit a truncated triangle 
side shape with typical width of 22.5 ± 0.12 nm and homogeneous height, of 2.4 nm (figure 

inb). 

In the present paper we show that the structure factor of embedded nanostructures can be 
directly determined by means of grazing incidence anomalous x-ray diffraction, allowing to 
recover the average height and strain of the QSs, to determine their composition and check 
the As/P exchange. As a matter of fact, strain of buried nanostructures is not directly 
related to composition (Vegards's law), but it also depends on size, morphology, and cap 
layer thickness. Then, tuning the x-ray energy near an absorption edge of atoms that 
belong to the nanostructures is a way to modify their scattering power and to enhance 
the chemical sensitivity of diffraction. Previous studies have shown the interest of using 
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Figure 1: Transmission Electron Microscopy images of InAs QSs embedded in InP (a) plane view, 
(b) cross section ; one can clearly see the truncated triangle shape 



anomalous diffraction to determine strain and composition in InAs QWrs grown on InP 
or large Ge dots grown on Si [3, We report on a general method that takes advantage of 
the full capability of anomalous diffraction and can be applied to the very interesting and 
challenging case of small size embedded nanostructures, the x-ray scattering yield of which 
is overwhelmed by the dominant matrix contribution, whatever the momentum transfer is. 

Grazing incidence anomalous diffraction intensity, can be written in the frame of the 
Distorted Wave Born Approximation Q, [l]| : 



I{Q,E)^\\T,{ai,E)\\' Ft{Q) [{cos{ipT - V a) + (3 f'As) + {^H^t - ^ a) + f'L) \ (1) 
where Ft is a complex structure factor of phase ipT which includes the overall contribution 
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of non anomalous atoms and the Thomson scattering of all anomalous atoms, (3 — \\p^\\ 
where Fa is a complex structure factor of phase if a which includes the Thomson scattering 
of all anomalous atoms (i.e. As atoms). ||rj(aj, .E) ||^ is the incidence transmittivity, and aj 
is the incidence angle that is close to the critical angle ac- In the following, the exit angle is 
much larger than the critical angle, then the exit transmittivity \\Tf{af,E)\\^ = 1 and the 
scattering length lOOnm) is much larger than the epilayer thickness. The As scattering 
factor writes fAs = Ias + Ias + '^fls^ where /^^ and /^^ are the real and imaginary resonant 
scattering of As atoms. Equation [T] shows that diffraction measurements at various energies 
(at least 3) at the As K-edge, allow to recover a quantity that is proportional to the Ft 
modulus, the /3 ratio and the phase difference A^p = lpt — Va- The knowledge of ||F4||and 



A93 readily gives information about the average size and strain of the nanostructures and 
chemical mixing at interface. 

Grazing incidence Anomalous diffraction at the As K-edge (11.867 KeV) was performed 
at the French Collaborative Research Group beamlines BM32 and BM2 at the European 
Synchrotron Radiation Facility, using a Si(lll) double crystal monochromator to select 
energies and mirrors to reject harmonics. Grazing incidence was used to minimize the 
substrate contribution. We recorded the scattering intensity in the vicinity of the (442) 
substrate reflection, at several energies across the As K-edge. We chose a weak reflection 
(h+k+l=4n+2) for which the anomalous diffraction contrast is maximum; it lies in the 
mirror plane defined by the [110] and [001] directions to benefit from symmetry in the 
reciprocal space. Figure El shows the diffraction intensity map recorded at 11.840 keV and 
at a grazing incidence angle equal to the critical angle (oj = ttc = 0.2°) to maximize the 
diffraction intensity. For that reflection the exit angle «/• is about 20°. A large amount of 
information can be drawn from such a map : the spreading of scattering in the [110] direction 
is due to both short range correlation and lattice strain in the corresponding direction in 
the real space ([110]) whereas, in the [001] direction it is due to the sharp strain evolution. 
Figure [21 shows correlation satellites on both sides of the (442) (SI and S2 ) their positions 
relative to the substrate peak give a qualitative estimation of the mean correlation distance 
between the wires. We find a value of 20. Trim, that is in agreement with the TEM. 

The markable result in this map, is the clear asymmetry of the SI and S2 satellites 
; SI exhibits a broad and twin feature along the [001] direction, with a sharp splitting 
at 1=2, whereas S2 is a single feature centered approximately at 1=2. SI corresponds to 
regions stretched along the [110] direction (tensile strain), i.e. to InAs sticks and InP regions 
located below and above the sticks. 82, instead, corresponds to InP which is compressed 
in between of the sticks. Indeed, anomalous diffraction measurements at the As K-edge 
show no significant intensity variations of satellite S2 as a function of the energy, confirming 
that the InP contribution is the dominant one. Furthermore, these measurements also show 
firstly, that the transmittivity corrections can be neglected, secondly, that the wetting layer, 
if it exists, must be very thin, at most one ML. Cutoff fringes are clearly observed in the / 
direction around the substrate Bragg peak position and give the total thickness of both the 
InAs QSs and the InP capping, that is equal to 11 ± O.Snm. 

We performed Z-scans {h=k=3.98) accross the satellite SI, at nine different energies close 
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Figure 2: Experimental diffraction map around the weak (442) reflection recorded at 11.840 KeV 
in grazing incidence geometry {oi = Oc = 0.2°, aj ~ 20°). SI and S2 are correlation satellites due 
to the stick short range periodicity. 
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Figure 3: (a) Experimental Ft (solid line) and Fa (ooo) modulus as a function of the reciprocal 
lattice unit I at /i=A;=3.98 (across the satellite SI). Fa is the structure factor that corresponds to all 
anomalous atoms, i.e. As atoms only. Also shown are best simulation curves (dashed line) obtained 
with a EDM model made of pure InAs QSs with a truncated triangle side profile. Inset shows 
the experimental (solid line) and simulated (dashed line) square root of the diffraction intensity at 
11.867keV ; at this energy, anomalous diffraction is maximized, (b) GI- DAES spectrum recorded, 
at the As K-edge, at the maximum of the Fa profile (Z=1.9) and the best fit curve obtained with 
pure InAs structure (solid line). 



to the As K-edge including those corresponding to the minimum of /^^ and the white line 
of f'^g. A weak fluorescence background was substracted and the data were normalized to l- 
scans across the satellite S2, measured at the same energies. Then Fa{Q), Ft{Q) and ^ip{Q) 
were extracted by fltting equation [T] to the energy dependent experimental intensities. 

Figure Et) shows the experimental modulus of Fa{Q), Ft{Q). The QSs height average 
value {{H)) can be estimated from the Full Width at Half Maximum (FWHM) of Fa , 

(H) = n = 2.54nm, that is very near of the value measured with TEM. In 

order to determine the QSrs composition and the local strain accomodation, Grazing In- 
cidence Diffraction Anomalous Fine Structure spectrum was measured at the maximum of 
Fa (h=k=3.98 and 1=1.9) at the As K-edge. Figure [Hb) shows the experimental DAFS 
spectrum and the simulation calculated with the crystallographic structure of InAsi-a;?! 
and the anomalous scattering factors /^^ and /^^of bulk InAs. A scale factor, the detector 
efficiency as a function of the energy and the As occupation factor (1 — x) were refined. The 
best fit curve, shown on figure Eb) corresponds to a value of (1 — x) equal to 1.03, i.e. the 
QSs composition is pure InAs. GIDAFS oscillations, in the energy range above the edge, 
will give also direct information on the local composition and on strain accomodation inside 
the sticks. Detailed analysis of these oscillations will be reported elsewhere. 

Finite Difference Method simulations were performed to map the strain produced by 
InAs QSs embedded in InP and compare its Fourier Transform to experimental diffraction 
intensity maps. We assumed a coherent growth of InAs on InP (supported by the lack of 
dislocations seen in TEM micrography). We performed Finite Difference Method calcula- 
tions in a 2D periodic frame, in the plane determined by the [110] and [001] crystal 
axis, i.e. the wires have infinite length along [110] and a periodic structure along the [110]. 
The QSs length is finite and rather short (50-200nm), meaning 2D calculations slightly over- 
estimate the strain, the relaxation along [110] being not permitted. We chose an FDM cell 
with a size of 4^ X f . Then, atoms are placed inside the crystallographic cell according to 
Cubic Face Centered Zinc-blende structure, taking into account elastic deformation. Start- 
ing from atomic positions, diffraction intensity is then calculated using the Distorted Wave 
Born Approximation [l^. No structural disorder has been taken into account. The P and 
In anomalous scattering factors were obtained from theoretical values, whereas experimental 
values for As in bulk InAs were used. 

Figures SK) andlUb), show the deformation maps with respect to InP lattice parameter, 



o 




(a) (b) 

Figure 4: FDM simulations of (a) £xx = "^'^aj"^'p^ and (b) £22 = '^'^"'^^"'^ , where aij and Cij are 
the parameters of cell (i,j). The strain is relative to bulk InP cell parameter, it is not the relative 
displacement. The Ezz cleary show the morphology of In As QWrs. Solid and dashed lines represent 
positive and negative contour strain, respectively, starting at —5.10"^ and spaced by .001. Ten 
additional Ezz contours are shown inside the wire starting at 0.06. The solid line on top represents 
the sample surface. 

obtained with a model of embedded InAs QWrs as deduced from TEM images (fig. [TJd), 
that well reproduces the diffraction data. The QWrs side shape is a truncated triangle as 
shown by TEM. The total thickness of the QWrs and the capping is equal to llnm, as 
reported above. Figure HJd) shows a map of deformations along the [001] direction {Szz), 
which emphasizes the morphology of the wires. Indeed the wires are compressed along [110] 
(with partial relaxation shown in fig. E^), thus expanded along [001]. This explains the sharp 
contrast of lattice strain shown in figure IHb) {Szz)- Figure S^) shows the relative strain Exx 
along [110], here it is impossible to disentangle InAs from InP regions under and on the top 
of the wires, since the in-plane lattice strain is continuous at the interface. On figure |Hti) 
are shown the best calculated curves of Fa and Ft, obtained by optimizing the height and 
width of the wires, as well as As/P intermixing at the InAs/InP interface. The splitting of 
Ft is well reproduced by the simulation and is due to a scattering phase shift between InP 



cells above and below the wires. The overall wire thickness and width are found to be 9 MLs 
and 22nm, respectively. The simulated ratio P = j-o^'^^|,|| at the maximum of is equal 
to 0.113, i.e close to the experimental value of 0.124 deduced from the GI-DAFS spectrum 
(figure Eb). Note that FDM simulation well reproduces the relative positions of Fa and Ft, 
with a strain Szz of about 6.3% in the inner part of the wire. For comparison, the strain 
of a pseudomorphic InAs thin film grown on InP, as foreseen by the elastic theory, is 6.7%. 
The e^z = ~ ^^"1^ value, deduced directly from the reciprocal lattice position Ipj, 

of the maximum of Fa (fig. Et) is equal to 6.1±0.25%, i.e., within the uncertainty, equal 
to the value simulated with FDM. Regarding composition, our FDM simulations show that 
the QSs inner part is pure InAs. The experimental curves {Fa and Ft) are compatible 
with a weak As/P intermixing at the InP interface, that would spread over one ML. Indeed, 
improving the signal-to-noise ratio together with a 2D data treatment, would allow to fully 
exploit the technique sensitivity to map the composition at the ML scale. 

In conclusion, we have shown that anomalous diffraction can be used to extract the 
structure factor of small size InAs nanostructures embedded in InP matrix. This study is a 
first step towards a 2D and 3D analysis of anomalous diffraction maps and Grazing Incidence 
DAFS data to recover the strain, size, shape and composition of embedded nanostructures. 
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